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An anomalous diffusion model for ion channel gating is put forward. This scheme is able to de- 
scribe non-exponential, power-law like distributions of residence time intervals in several types of 
ion channels. Our method presents a generalization of the discrete diffusion model by Millhauser, 
Salpeter and Oswald [Proc. Natl. Acad. Sci. USA 85, 1503 (1988)] to the case of a continuous, 
anomalous slow conformational diffusion. The corresponding generalization is derived from a con- 
tinuous time random walk composed of nearest neighbor jumps which in the scaling limit results 
in a fractional diffusion equation. The studied model contains three parameters only: the mean 
residence time, a characteristic time of conformational diffusion, and the index of subdiffusion. A 
tractable analytical expression for the characteristic function of the residence time distribution is 
obtained. In the limiting case of normal diffusion, our prior findings [Proc. Natl. Acad. Sci. USA 
99, 3552 (2002)] are reproduced. Depending on the chosen parameters, the fractional diffusion model 
exhibits a very rich behavior of the residence time distribution with different characteristic time- 
regimes. Moreover, the corresponding autocorrelation function of conductance fluctuations displays 
nontrivial features. Our theoretical model is in good agreement with experimental data for large 
conductance potassium ion channels. 

PACS numbers: 05.40.-a, 87.10.+e, 87.15.He, 87.16.Uv 



I. INTRODUCTION 

Ion channels are complex membrane proteins which provide ion-conducting, nanoscale pores in the biological mem- 
branes [lj. These proteins undergo spontaneous conformational dynamics resulting in stochastic intermittent events 
of opening and closing the pore - the so-called gating dynamics. It can be described by following kinetic scheme: 

C^O. (1) 

k c 

As it stands, this scheme describes Markovian stochastic transitions between the closed state (C) and the state open 
(O) of an ion channel which can fully be characterized by the opening rate, k a , and the closing rate, k c . From a 
trajectory description of the observed two-state gating process, these transitions can be characterized by the residence 
time distributions (RTDs) of open and closed time intervals, iPo(t) — fc c exp(— k c r) and ip c (T~) = k exp(— k r), 
respectively. 

The invention of patch clamp technique Q marked the beginning of a new era: detailed experimental studies of the 
statistics of such stochastic trajectory realizations have been rendered possible. These experimental investigations, 
however, also reveal the fact that the distributions of the residence time intervals are typically not exponential. This 
in turn implies that the corresponding observed two-state dynamics of current fluctuations is not Markovian. Any 
such non-exponential distribution can however approximately be fitted by a sum of (sometimes many) exponentials, 
e.g. 

N 

^c{t) = ^2 exp(-Ajt), (2) 

i=l 

with weights Cj obeying X^iLi c » = 1- The rationale behind this fitting procedure is the assumption that the corre- 
sponding state consists of N discrete substates, separated by potential barriers. This method constitutes the working 
tool for the majority of molecular physiologists in interpreting their experimental data within a discrete Markovian 
scheme consisting of many (sub)states Q. The addition of new states (or new configurational dimensions in the 
continuous case) is a well known formal method to unravel a low-dimensional non-Markovian stochastic dynamics 
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via its embedding into a Markovian dynamics of higher dimensionality. The problem with such a methodology is, 
however, that the number of substates needed to fit the experimental data can depend on the experimental conditions. 
For example, the experimental gating dynamics of a Shaker potassium channel has been successfully described by 
a sequential 8-state Markovian scheme with 7 closed states for a fixed value of temperature about T = 20 °C [j]. 
However, to describe the experimental data over a small extended temperature regime between 10 — 20 °C already 
necessitates to add three additional closed substates p|. 

For several types of ion channels the RTDs can alternatively be fitted in terms of non-exponential distributions such 
as a stretched exponential i^{t) oc — exp[— {vr) a \ 0, or by a power law V'( T ) °^ with a few parameters 
only. The case of a power law coefficient near (3 = 3/2 can be described with normal conformational diffusion over 
many degenerate substates H US [H d. Such diffusion models HUHEJHEiElHIll P resent an 
alternative to the standard discrete state Markovian modeling which can be considered also as a complementary one 
|l6j. It should be noted, however, that exponents (3 different from a normal diffusion behavior have been detected 
experimentally as well [I7l . ll8l.ll9j . Therefore, it is of prominent importance to generalize the normal diffusion model 
to the case of anomalous diffusion: this objective is at the heart of the following discussion. 



II. MODELING ANOMALOUS ION CHANNEL GATING 



Ion channels are protein complexes with an intrinsically hyperdimensional conformational space. Such macromolec- 
ular complexes can possess several macroconformations corresponding to different functional states of the ion channel. 
These macroconformations are, however, not static, but rather present dynamical structures featuring multiple con- 
formational substates and undergoing a corresponding intraconformational dynamics which can be of the diffusion 
type, as it is well studied for one of the simplest proteins - myoglobin [2(], [^J [^ . It particular, for myoglobin it has 
been shown experimentally that the diffusive motions become increasingly important above some critical temperature 
and can serve as lubricants for macroconformational changes, see e.g. in Ref. |22j . There is no convincing reason to 
believe that ion channels, which are structurally much more complex proteins, are very different from myoglobin in 
this respect. Therefore, the conformational diffusion should be present and can be important for the gating dynamics 
of ion channels as well. Diffusion models of ion channel gating put emphasis on this fundamental feature of protein 
dynamics. In the case of ion channels with one open and one closed macroconformations the whole conformational 
space can be divided into the two hyperdimensional subspaces corresponding to the open macroconformation and to 
the closed macroconformation, respectively. These subspaces present generally some complex manifolds consisting of 
domains of attraction in the conformational space which are separated by potential barriers. Such domains of spatial 
localization can be associated with discrete states in the discrete state diffusion modeling of the gating dynamics. 
They can possess but a complex intrinsic structure reflecting the intrinsic multidimcnsionality of the problem. For 
this reason, such states can be trapping states with a non-exponential distribution of the residence times spent in the 
corresponding domains of attraction. On the contrary, the Markovian assumption within a discrete state modeling in- 
volves the assumption on strictly exponential distribution of the residence times in such discrete states. This presents 
at best an approximation only. It is at the heart of modeling the gating dynamics with appropriate Markovian kinetic 
schemes. Moreover, some parts of conformational potential landscape can lack the presence of pronounced barriers 
and some domains of attraction can be flat providing valleys in the corresponding potential landscape. This provides 
an ideal starting point for a diffusion modeling of the gating dynamics. 

The simplest diffusion model El assumes the existence of such a one-dimensional valley wherein the diffusion is 
modeled by sequential transitions am ong a large number of discrete substates. Alternatively a continuous diffusion 
modeling can be applied [UJ, |ll|, Il4. Il5| . Note also that multiple pathways into the open state can be accounted 
for within a one-dimensional reaction coordinate modeling by allowing for transitions with non-nearest neighbor 
jumps. Such complexity will not be addressed within our sequential diffusion model. At one edge of this valley the 
transition into the open state can occur via the single route El E3 ■ Such a modeling is appropriate when the 
correlations between the subsequent closed and/or open residence time intervals are absent. The multiple transitions 
could also induce some correlation between the durations of time intervals. Such effects will be neglected in a zero 
order approximation (renewal assumption for the observed two-state conductance fluctuations). The existence of 
time-correlations between the durations of time intervals (non-renewal gating dynamics) would indicate the presence 
of more than one route between the manifolds of open and closed states (see, e.g., in Ref. Q, PP- 440-441). Such 
correlations can also be accounted for in the diffusion models 0, The incorporation of such correlations is, 

however, beyond the scope of our present work which presents an anomalous diffusion generalization of the Markovian 
modeling used in Ref. [71 l9l Hoi m| . 

Let us start from a continuous time random walk (CTRW) [13, 0, HfJ |2(| |2]j generalization of the discrete state 
diffusion model by Millhauser et al. 0, see Fig. 1. It is assumed that the manifold of closed substates consists of N 
states; namely, the states from j = 2 to j = N — 1 ("diffusion states") are identical and characterized by identical 
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FIG. 1: Sketch of the CTRW generalization of the (discrete) diffusion model of ion channel gating. 

residence time distributions tpj(r) = ip(r), i.e. the channel stays in the corresponding state j for a random time 
interval r distributed in accordance with ip(r), and performs at the end of every time interval a jump j — ► j ± 1 
with probability Pj+i,j = Pj-i.j = 1/2 cither to the left, or to the right neighboring state, respectively. If the RTD 
is exponential, i.e. iJj(t) — 2«;exp(— 2kt), then the standard Markovian rate description with rate k is recovered. 
The boundary state j — N possesses a different RTD iPn(t) which in the Markovian case reads i/jn(t) = Kexp(— kt) 
(transitions occur always with the probability one, pn-i,n = 1, to the state j = N — 1). Furthermore, from the 
state j = 1 the channel can undergo a transition into its open state j = with the rate 7, or make transition into 
the manifold of conformational diffusion substates with the rate n. For this state, the corresponding RTD reads 
ipi (t) = (7 + «)exp[— (7 + k)t] and the transition probabilities are poi — 7/(7 + k) and P21 = k/(j + k). The 
dynamics of state occupancies Pj(t) is described by the generalized master equation (GME) due to Kenkre, Montroll 
and Shlesinger j2^, |2£j and its generalization [3(3, |ijnj • The corresponding dynamics reads (with an initial preparation 
in some state jo, Pj o {0) = 1): 



Pj(t) = / K(t-t')\p j -i(t')+p j+1 (t')-2p j (t')]dt\ j = 3,N-2, (3) 
Jo 

p N {t) = f ^(t-t')Piv-iW- / K N (t-t')p N (t')dt', (4) 

Jo Jo 

PN-i(t) = f K(t-t')[p N _ 2 (t')-2p N _ 1 (t , )}dt' + f K N (t-t') PN {t')dt', (5) 

Jo Jo 

P2 (t) = [ K(t-t')\p 3 (t')-2 P2 (t')}dt' + K Pl (t), (6) 







t 



pi(t) = / K(t-t')p 2 (t')dt'-(K + j)pi(t) + k c p (t), (7) 
Jo 

po(t) = ~k cPo (t)+ 1Pl (t), (8) 
with the kernels K(t) and Kpj(t) defined through their Laplace-transforms 

K{s) = L s ^ 



2 l-^(s)' 
Zr 1 \ s ^n{s) 

Kn{s) = r — , (9) 

1 - tpN(S) 

where ip(s) and iPn(s) are the Laplace-transforms of i/>( T ) an( i ^n{t), respectively. These kernels can be derived due 
to the following simple consideration. Let us prepare the diffusing particle at t = in the state j = 2...N — 1 and let 
it make transitions to the neighboring states without return. Then, the survival probability <£>(£) — ip(r)dT in the 
state j (which should be identified under such conditions with Pj(t)) is governed obviously by the following equation 
$(i) = -2 f* K(t - t')$(t')dt' with $(0) = 1. It can be easily solved by the Laplace-transform method to yield the 
first equation in Eq. J!2J) by taking into account ip(s) = 1 — s$(s). The second kernel K^(t) is obtained along the 
same lines. 

The RTD of the open state is readily obtained, i.e. ipo(i~) = fc c exp(— k c r). In order to calculate the RTD of the 
set of closed states one starts out from pi(0) = 1 (the channel has been just closed) to obtain the survival probability 
^c(i) = J2j=iPj(t) with the boundary condition that the state j = is absorbing. This latter condition is realized 
by setting formally k c — > 0. The corresponding RTD then follows from ^> c (t) = — d$ c (T) /dr. The total population of 
the closed state p c {t) = YljLiPj(t) obeys (not allowing for the backward transition, k c — > 0): 

Pc(t) = -lPi(t). (10) 
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This must be used as the proper boundary condition (it yields a radiation boundary via continuity equation in the 
scaling limit, see below) to calculate ^> c (t) |33ll34||. 

Next let us make the Ansatz that t/>(r) = --£:E a (-(2KT) a ), where E a (z) := z n /T(an+l) is the Mittag-Leffler 
function. It is defined via a generalization of the Taylor series expansion of the exponential function, E\{z) = cxp(z) 
and r(z) is the standard Gamma-function. In other words, the corresponding survival probability in the state just 
occupied gjj, $(r) = f°° ip(t)dt is given by $(r) = E a (— (2kt) q ). This corresponds to the Cole-Cole relaxation 
law, i.e. the relaxation law which yields the Cole-Cole susceptibility Isfij. The Laplace-transform of $(t) reads 
$(s) = s a - 1 /[a a + (2k) q ] [HI43 and by use of the relation, $(s) = 1 - s$(s), one obtains ^(s) = (2«;) a /[s" + (2k)"]. 
This particular choice of RTD interpolates between the initial stretched exponential (Weibull) distribution |4jj (which 
corresponds to the Kohlrausch relaxation law; "stretched exponential" refers here to the survival probability $(r), 
with ip{ T ) l/r 1_Q at r — > 0) and the asymptotic long time power law ip( T ) <x 1/t 1+q . It yields anomalously slow 
diffusion, {5x 2 (t)) oc t Q 42]. For example, such an anomalous diffusion is measured experimentally, along with the 
RTD in trapping domains exhibiting the corresponding power law, for colloidal particles in cytoskeleton actin networks 
of biological cells [3jj and in living cells |38j . These latter experimental results offer a clear clue for understanding 
the results on virus diffusion in infected cells |39| - that is an observed anomalously slow diffusion in actin networks 
combined with active directional trafficking of viruses by molecular motor proteins 39] . For the discussed form of 
RTD, K(s) — (2k) q s 1_q /2 and the fractional master equation follows exactly from Eq. J3J as a particular case of the 
GME. It reads explicitly, 

pj(t) = i(2«)« ^~ a [pj-i(*') +Pj+iV) - 2 Pj(t% (11) 

where oD t (...) = fXajlk Jo (t- t '\^ r s 1S ^ ne integro-differential operator of fractional derivative introduced by Rie- 
mann and Liouville, see Refs. [4ll l42t 1431 l46j for reviews and further references. In the case of a two-state dynamics a 
similar fractional master equation was obtained in 0, . Note that the fractional master equation presents in 
fact a conventional generalized master equation being nonlocal in time. The introduction of a fractional time deriva- 
tive in the generalized master equation of CTRW is nothing but a shorthand notation which corresponds to a specific 
choice of the RTD. The importance of this equation lies in the fact that it can serve as a useful mathematical tool 
to model anomalously slow diffusion. The physical origin of this diffusion can be attributed to very broad residence 
time distributions on the sites of particle localization with diverging mean residence time [23L l25l |2(j . In practice 
this implies that the corresponding mean residence time is exceedingly large as compared with the characteristic time 
scale of anomalous diffusion in the given domain of a finite size. As a consequence, the approximation with an infinite 
mean residence time becomes physically justified. 

Likewise, with ^jv(t) = — -^E a {~{nT) a ), Eq. takes on the form 

p N (t) = ~(2n) a D t ~ a p N -i(t') - K a D~ a PN (t r ). (12) 

The remaining equations JSJ), and J7J involving the memory kernel can readily be rewritten in a similar form 
upon use of the notation of the fractional time derivative. 

It is worth noting that Eq. (1111) and Eq. (|12fl can equivalently be brought into a form with the fractional Caputo 

derivative £)"(...) := r(i-a) Jo E3 m the left hand side. Namely, Eq. (|11|) becomes 

D>At') = \{^T\Pj-x{t) +Pj+i(t) - 2pi(t)] • (13) 

Such a fractional master equation has been first derived from CTRW in Ref. j4^- This form is, however, clearly not 
appropriate for a further use in the studied case, where the normal rate transitions are also present. The reason is 
that the Eqs. 10, would get a form mixing the fractional Caputo derivative and the fractional Riemann-Liouville 
integral in left and right hand sides of the corresponding equations. This fact establishes a clear advantage for the 
use of the Riemann-Liouville fractional derivative in a case where normal rate processes are also present. 

A. Scaling limit to a fractional diffusion equation 

Let us perform next a continuous limit: namely, assuming the distance Aa: between neighboring sites we introduce 
the conformational coordinate x := —jAx which models the manifold of closed diffusion substates. The following 
continuous limit is assumed: Let Ace — > 0, N — > 00, k — > 00 whereas keeping K a :— ^(2K) a (Ax) 2 and the diffusion 
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"length" L := NAx constant. By use of the expansion: p j±1 (t)/Ax := P(-\j ± l]Ax,t) « P{x,t) =F aP ^ :t) Ax 



- 3 p ( a ' t ) (Ax) 2 , in Q we arrive at the following fractional diffusion equation in continuous state space, i.e., 



9P(x,t) fl-a d*P(x,t') 

—Q^=K a0 D t ^ 2 , (14) 

where K a is the diffusion constant of anomalously slow diffusion. The fractional diffusion equation l|14|l presents a 
non-Markovian diffusion equation |2?j with a particular choice of an integral kernel that accounts for the residence 
time distributions following the Cole-Cole relaxation law |35Ll36| in a corresponding formulation in terms of a discrete 
CTRW. The fractional diffusion equation (|14|l assumes the form of a continuity equation, 

dP(x,t) dJ(x,t) 

~^r~ = — dx~ (15) 

where the probability flux J(x, t) becomes modified due to the fractional time-derivative, 

~ l— a dPix, t') , „. 

J(x,t) = -K a0 D t — ■ ( 16 ) 



Our derivation of the fractional diffusion equation from the CTRW complements previous studies |42|, |4JJ ; it is rather 
simple and does not require jumps with a variable step length beyond nearest neighbors. For this very reason, no 
overflights of the boundaries occur that are possible otherwise. This observation is of crucial importance in determining 
the physically correct boundary conditions 0, ■ This also means that the boundaries are strictly local in space. 
The given derivation removes any open query about the space-locality of boundary conditions for the fractional 
diffusion equation. After the integration of the continuity equation from x = —L to x = 0, one deduces that the 
decrease of the total probability of the closed state manifold (the survival probability), p c {t) = J° L P(x, t)dx occurs 
due to the probability flux on the boundary. Accordingly, we will replace the original discrete master equation in 
space by its corresponding fractional diffusion equation with the following boundary conditions as they emerge from 
the original problem. The boundary x = —L is a reflecting one, obeying : 

J(-L,t) = 0. (17) 

The boundary at x — is radiative. Setting p%(t) w AxP(0,t) = LP(0,t)/N and using Eq. (|10fl one finds 

J(0,t)=jLP(0,t)/N. (18) 

We additionally use here the specific scaling limit: 7 — > 00, N — > 00 with (r c ) = N/j being held fixed. This quantity 
possesses the meaning of the mean residence time in the closed state manifold, see below. The radiation boundary 
acquires then the explicit form 

J(Q,i) = LP(0,t)/(T c ). (19) 

Our fractional diffusion modeling for the RTD of closed time intervals thus has three parameters only: (i) the mean 
residence time (r c ), (ii) the characteristic time of conformational diffusion, i.e. td := (L 2 /K a )^ a and (iii) the power 
law index of anomalous diffusion a. In the case of normal diffusion, i.e. a = 1, the special scaling limit used here can 
also be justified from a Kramers approach [5l| to the gating problem 0,0]. Whereas the solution of the Kramers 
approach in Ref. 11] in addition yields also an analytical expression for (r c ), which reproduces the experimental 
crossover behavior from an exponential-to-linear voltage dependence due to Hodgkin and Huxley |5^ | , the model here 
treats (r c ) as one of the phenomenological parameters. It is also worthwhile to remark that the boundary condition 
does not contain the index of anomalous diffusion a. Formally it remains the same as for normal diffusion. The flux 
expression (I16fl is, however, not local in time. Moreover, the r.h.s. of (|19|) does not contain the fractional derivative 
in time. This feature is in accord with the original discrete model where the last, final transition into the open state 
is given by an ordinary rate transition. The analogy with the Kramers model of Ref. |ll| is that the diffusion in 
the domain of voltage-dependent states (cf. [ITIIT^ ]). which becomes a thin boundary layer in the considered scaling 
limit, remains normal. This justifies well the use of boundary condition (|19fl in our fractional diffusion model which 
now is completely formulated. 



B. Characteristics of the residence time distribution 



To obtain the distribution of closed residence times iPc{t) one needs to solve first the fractional diffusion equation 
(f 141) with the boundary conditions (|17fl and (|19|l and the initial condition P(x,0) = 5(x — xq) with xq — > 0_. 
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The survival probability $ c (t) follows as the integral of the solution over the spatial variable and subsequently the 
corresponding RTD follows as the negative time derivative of the survival probability. This task has been achieved 
by use of the Laplace-transform method. The details of the derivation are outlined in the Appendix. The final result 
for the Laplace-transformed RTD of closed time intervals then reads: 



where an auxiliary function 



Ms) = - I r , (20) 

1 + s(T c )g a (sT D ) 



, . tanhfz"/ 2 l , s 

M*) = L 1 (21) 



has been introduced. For a = 1, this result reduces to one for normal diffusion in Ref. [Tllll2j |. Moreover, since g a {z) = 
1 + o(z) for small z, one can readily see from Eq. I|20[l that (r c ) indeed has the meaning of a mean residence time, 

(t c ) := J °° t^P c (t)cIt — — d ^ d ^ | s — o- Note also that the second moment of RTD diverges, (t 2 ) := T 2 tfj c (T)dT — ► oo 
for all a < 1 (anomalous diffusion). Furthermore, if td = 0, then the closed time distribution becomes strictly 
exponential, i.e. ip c (r) = cxp(— t / (t c ))/ (t c ) and the simplest two-state Markovian model of the ion channel gating is 
reproduced with the opening rate k Q = (r c } _1 . In general, the expression l|20|) cannot be inverted to the time domain 
exactly; its different characteristic regimes, however, can be discussed analytically. 

In proceeding, let us consider first the limit of a large conformational diffusion time td 3> (t c ). Then (by use of 
the large- z asymptotic behavior of g(z) ~ z~ Q / 2 ), we have 

^ " 1 + (22) 

// ,\l/(l-a/2) 

with To := td[ ■ The inversion of Eq. (|22ll is given as ip c (r) — ~d^ c (r) / 'dr in terms of the survival 

probability 

$ c (r)=B 1 _ a/3 [-(-J J, (23) 

which is expressed through the Mittag-LefHer function E a (z) and corresponds to the Cole-Cole relaxation law [85ll3^ |. 
Because Ei/ 2 {— z 1 ^ 2 ) = e z erfc(z 1//2 ) [42| . where erfc(z) is the complementary error function, the solution of the normal 
diffusion problem in Ref. [10| for the initial and intermediate time evolution regimes is reproduced from Eq. 123fl . 
For t<t , Eq. (|2*3l behaves as a stretched exponential |4^|. 

This dependence 124|) corresponds in the language of time-dependent rates k (r) := — -4- ln[$ c (r)] (used in the renewal 
theory) @, E3 to fe (r) oc T~ a l 2 . Such a time-dependent rate of recovery from inactivation has been measured 
with a = 1 for a sodium ion channel in Ref. [Hij . In the limit t — » 0, Eq. (|24|l yields a power law for the RTD, 

Mr) oct-"/ 2 . (25) 

Furthermore, for intermediate times r <C t Ctjj, Eq. I|23|) yields another power law, reading 

Mr) oc (26) 

with [3 = 2 — a/2. For a — 1 such an intermediate power law with the slope (3 — 3/2 has been measured for a 
potassium ion channel in Note that this particular power law exponent reflects normal diffusion. Consequently, 
the origin of the intermediate power law is principally not due to the anomalous diffusion behavior. Our theory 
smoothly reproduces the intermediate power law associated with normal diffusion in the limit a — ► 1. 

Other power law features were also measured in experiments, for example, for a gramicidin channel with the slope 
(3 w 1.7 |l7j . This corresponds to an intermediate power law in Eq. I)26ll with a « 0.6. Moreover, power law exponents 
with (3 > 2 are also measured in experiments |l<Sl analyzed in Ref. |19| from a pure phenomcnological perspective 
without clarifying a tentative mechanism. Our model can as well capture such anomalous power laws which cannot 
be explained within the intermediate power law asymptotics (|26ll . 
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Indeed, let us assume that Td is sufficiently small and to consider the asymptotic behavior r — > oo. The corre- 
sponding asymptotics can be deduced from the behavior of V> c (s) at small s. For s — > 0, Eq. (|20|) yields, 

^ c ( S )«l-s(r c )[l-( S T D ) Q /3]. (27) 

From this, by way of $ c (s) = (1 - V> c (s))/s » (r c )(l - (st^)"^) for s ->■ 0, it follows that $ c (r) oc 1/t 1+q for 
t — > oo. This renders then a power law (|26|) with (3 = 2 + a for large r. This asymptotic power law with (3 > 2 
is a manifestation of the anomalously slow conformational diffusion in a space domain of finite size. It replaces an 
exponential asymptotic behavior of ip c (t) for r > td in the case of normal diffusion 0,EI0- 



III. APPLICATION TO GATING DYNAMICS OF A LOCUST POTASSIUM ION CHANNEL 

Our fractional diffusion model can be used to describe the rather complex gating behavior observed for a locust 
potassium ion channel [l9j . This ion channel exhibits experimentally a Pareto law in its gating kinetics, ^> c (t) = 
a/ (& + r) /3 with (3 ~ 2.24 ±0.06. The corresponding autocorrelation function, however, seems to exhibit three different 
interchanging power laws |l9j |. These features are compatible with our model. Within a two-state reduction we are 
dea ling with an alternating renewal process [40l | . Its (Laplace-transformed) normalized autocorrelation function reads 



1 f 1 1 \ 1 



i - i> (s)) (l - 

k ^ = 1 -[7—\ + 7—\)-—/ \ — (28) 

s \(t„) {t c )J s 2 (l_^ ( s)V , c ( s )j 



For our case under study this yields 



fa(sT D ) + S(T C ) 



*(«) = - (gr . : \ , r . ( 29 ) 



s l + ^ + f a ( S r D ) + s(T c )' 



where f a (z) :— l/g a (z) — 1. Note that / a (0) = and for td — the inversion of Q29f) yields the Markovian result 
k{t) = exp[— (ho + k c )t]. Moreover, the analytical expression (|29|l allows one to study the asymptotics of k(t) at 
t — > oo. Namely, from / q (std) ~ (st£>)"/3 it follows that k(s) oc s Q_1 at s — > 0. By virtue of a Tauberian theorem 
[58| . this latter result readily yields, 

k(t) oc r a . (30) 

This power law feature agrees well with the experiment which shows asymptotically H30fl with a — 0.28 ±0.1. 
Furthermore, an intermediate asymptotics of k(r) can be obtained by studying the limit of very large td- Using the 
scaling s :— s(t c ) and the limit of very large values y := td/(t c ) 3> 1 such that s 1 is allowed for, whereas still 
sy 3> 1, Eq. (|29|l can formally be approximated by 

W2-1 

with r = r^ 1 (l + (r c )/(r )) 2 / a . The formal inversion of Eq. (2D yields k(t) = E a/2 [-{rt) a / 2 ]. This in turn yields 
an intermediate asymptotics k(r) oc r~ a / 2 within (r c ) < t < m. Indeed, the analysis of experimental data in 0] 
reveals such an intermediate asymptotics fc(r) oc T -°- 14 ± - 02 . The numerical inversion of k(s) in Fig. 2 displays three 
different power law regimes in qualitative agreement with the experimental data. Only one of these power laws - the 
long-time asymptotical one - seems, however, to present a true power law asymptotics. The intermediate power law 
in Fig. 2 does not agree numerically with the experimental one in Ref. [ljj- Nevertheless, the experimental data 
agree - surprisingly enough - with the intermediate asymptotics obtained above in the limit td — y oo. 

Furthermore, the numerical inversion of ipc{s) m Fig. 3 can be fitted by Pareto law with /3 « 2.24. The discrepancy 
between [3 — 2 w 0.24 and a = 0.28 is due to the experimental restrictions on the maximal time intervals measured. 
The actual power law asymptotics for r — ► oo in Fig. 3 is ^ c (t) oc t~ 1,28 . This long time asymptotic regime is not 
yet attained in Fig. 3 which instead closely agrees with Q c {i~) oc t -1 - 24 , see in Fig. 3, giving an apparent agreement 
with the experimental data. In view of our few elementary model assumptions, the agreement between theory and 
the experimental data |l8j analyzed in Ref. [l9| is striking indeed. 

Our fractional diffusion scheme is not expected to describe the experimental facts quantitatively in all details. In 
particular, it predicts that the low-frequency part of the spectral power S(f) of ion current fluctuations of the locust 
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FIG. 2: Normalized autocorrelation function of conductance fluctuations. Numerical inversion of Eq. 129^ is done with the 
(improved) Stehfest algorithm 60] for the following set of parameters: (r c ) = 0.84 msec, (r ) = 0.79 msec [19j and assumed 
td = 100 msec and a = 0.28. 
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FIG. 3: Survival probability of the closed state for the studied model. The set of parameters is the same as in Fig. 2. 



BK potassium ion channel corresponds to l// 7 noise with 7 = 1 — a ~ 0.72 |57j,|59j. The experiment |6l| indeed 
reveals l// 7 noise with 7 close but to unity, 7 « 1. The reason for this discrepancy is not resolved. The asymptotic 
behavior of the autocorrelation function in Ref. |19j and the behavior of the low- frequency part of the spectrum in Ref. 
[6l| are certainly at odds. A possible reason could be nonstationarity of the given current recordings. Nevertheless, 
the qualitative agreement, i.e., the principal theoretical prediction and the measurement of l// 7 noise, is comforting. 

Moreover, the durations of residence time intervals in open and closed states can be correlated. Such correlations can 
be induced by stochastic binding of calcium ions which regulate the gating dynamics of large conductance potassium 
ion channels. To account for such correlation effects, our model principally can be generalized in the same spirit like 
the original diffusion model has been generalized to include ligand binding effects |62|. This generalization is left but 
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for a future study. 



IV. DISCUSSION AND CONCLUSION 



The gating dynamics of protein ion channels in biological membranes is governed by a conformational dynamics on 
a very complex energy landscape with a huge number degrees of freedom. This multidimensional energy landscape can 
possess deep energy wells (as compared with the thermal energy ksT) which are separated by potential barriers. In 
addition, there can exist an underlying energy valley network connecting these wells which results in an energy quasi- 
degeneracy. The traditional discrete state approach to the gating dynamics pursued by the community of molecular 
physiologists presents an abstraction to this complexity: it has its focus on the fact of deep potential wells being 
separated by high energy barriers. The energy quasi-degeneracy of potential wells enters the theory as an entropic 
contribution to the corresponding free energies after reduction of the multidimensional reality to low-dimensional 
models (possessing a few discrete states only). This traditional approach has proven useful over the past years and 
it serves as a serviceable working tool for the analysis of the experimental data. This approach is, however, not able 
to capture the physical origin of such complexity features as the presence of power law distributions of the residence 
time intervals, the slow decay of the autocorrelations of fluctuations or the presence of l// 7 noise feature in the power 
spectrum of fluctuations of several ion channel types, to name but a few. Experiments have demonstrated 63] that 
the l// 7 noise is due to the conformational transitions among different conductance states. In particular, the ion 
current is free of l// 7 noise in a frozen conductance substate of the ion channel. Therefore, the l// 7 noise originates 
due to fluctuations among experimentally distinguishable substates |63| . These features reveal unambiguously the 
non-Markovian character of the observed "on-off" ion current fluctuations [6J, I64L Ififij . The diffusion models of ion 
channel gating 0, U E3, O, 0, 0, 0, El present another, complementary abstraction of the actual dynamics. 
These approaches attempt to capture the spatial structures of the potential minima and the associated dynamics, 
and/or the corrugated and hierarchical features of the real multidimensional conformation landscape after performing 
the reduction to a reaction coordinate picture |20|. It is physically likely that the ion channel protein can become 
temporarily trapped in some domains of its intrinsic conformational landscape from which it can escape by activated 
jumps among those states. Due to a complicated structure of such traps the corresponding residence time distribution 
can possess a divergent, - or from a practical point of view - a very large first moment. This in turn gives rise to an 
anomalous conformation diffusion within the chosen reduced reaction coordinate description. Our scheme in terms of 
a fractional, continuous diffusion model, being complemented with appropriate boundary conditions properly accounts 
for such complexity. 

As demonstrated theoretically and exemplified with the gating of a BK locust potassium ion channel our fractional 
diffusion theory presents a powerful approach to describe these various observed power law characteristics of the 
underlying gating dynamics. 
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APPENDIX A: SOLUTION OF THE BOUNDARY- VALUE PROBLEM 



The Laplace-transformed probability P(x, s) := J °° e st P(x,t)dt, Eq. l|14fl reads: 



sP{x, s) — S(x — xq) = K a s 



^ a d 2 P(x,s) 
dx 2 



(Al) 



with — L < xq < 0. Note that the limit xq — + 0- will be taken at the very end of calculation. The corresponding 
Laplace-transformed boundary conditions assume the form 



dP(x,s) 



dx 

^_ a dP(x,s) 
dx 



x— — L 



0. 



x=0 



L 

K a (r c ) 



P(0,s) 



(A2) 
(A3) 
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The challenge is thus the solution of the boundary- value problem (!Al|l - ilA3f) . Towards this goal we consider separately 
the solution in the domains, — L < x < x , 



Pi(x, s) = Ax exp \ J j^-xj + Bx exp { — J J^~ x ) - 



and xq < x < 0, 



P 2 (x, s) = A 2 exp (W j^xj 



+ B 2 exp - a — 



At x = xq, the solution is continuous, 



Px(x ,s) = P 2 (x ,s). 



(A4) 



(A5) 



(A6) 



The first derivative dP(x, s)/dx, however, experiences a jump. This can readily be seen upon integrating Eq. I jAlj l in 
an infinitesimally small neighborhood of x = xq. Thus, 



K a s l 



dP 2 (x,s) dPx(x,s) 



dx 



dx 



= -1 



(A7) 



The coefficients Ax. 2 and B 12 are determined by substitution of l|A4(l and (IA5fl into Eqs. (|A2|) . (|A3ft (|A6f) and 
(|A7I) . Thereby, the proposed objective is exactly solved. The integration of the solution (|A4I) from x = —L to 
x = (xq — > Q_) yields the Laplace-transformed survival probability <& c (s); the corresponding RTD l|20|) follows as 
4> c {s) = 1 - s$ c (s). 
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